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Abstract 

We study discrete curvatures computed from nets of curvature lines on a given smooth surface, 
and prove their uniform convergence to smooth principal curvatures. We provide explicit error 
bounds, with constants depending only on properties of the smooth limit surface and the shape 
regularity of the discrete net. 

1 Introduction 

The field of discrete differential geometry has brought to light intriguing discrete counterparts of clas- 
sical differential geometric concepts as well as efficient geometric algorithms (see, e.g., Bobenko et al., 
2008; Grinspun and Desbrun, 2006). One aspect of this theory is convergence: classical smooth no- 
tions should arise in the limit of refinement. Recently, several convergence results have been obtained 
for curvatures and differential operators defined on polyhedral surfaces. Roughly, one may distinguish 
three approaches: (i) polynomial surface approximation (see, e.g.. Meek and Walton, 2000; Cazals and 
Pouget, 2005), (ii) geometric measure theory (see, e.g., Fu, 1993; Cohen-Steiner and Morvan, 2006), 
and (iii) finite element analysis (see, e.g., Dziuk, 1988; Hildebrandt et al., 2006). Among these, (i) 
provides pointwise convergent curvatures for many, but not all, discrete meshes. In contrast, (ii) and 
(iii) consider generalizations of integrated, or total, curvatures and yield convergence in the sense of 
measures or appropriate Sobolev norms, respectively. 

Given the convergence of curvatures studied by approaches (ii) and (iii) in an integrated sense, it 
is natural to ask whether these curvatures can be shown to also converge in a pointwise manner. An 
affirmative answer can be obtained in some special cases, such as polyhedral surfaces with vertices on 
the unit 2-sphere (Xu, 2006). In general, however, the answer to this question is negative: it was ob- 
served in Xu et al. (2005) that for general irregular polyhedral surfaces, there exist no k-local definitions 
of discrete curvatures that are pointwise convergent. Here, by /:-locality we mean that the definition of 
curvatures associated with a vertex of a polyhedral surface only depends on the A;-star of p, i.e., those 
vertices that are connected to by a path of at most k edges. The concept of /:-locality is motivated by 
the smooth setting, where the definition of curvatures and differential operators only depends on local 
properties of the underlying Riemannian manifold. 

Uniform convergence from nets of curvature lines We provide an affirmative answer to the above 
question of pointwise convergence of curvatures for a special class of discrete meshes: discrete nets of 
curvature lines on a given smooth surface M that is immersed into Euclidean space (see Figure 1). 
To obtain approximations {ki , k2) of principal curvatures (a:i , K2) on M, we follow a three-step approach. 
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Figure 1 : A discrete net of curvature lines on an ellipsoid. 



We consider (a) local polyhedral approximations to nets of curvature lines, to which we apply (b) well- 
known 1-local integrated notions of discrete curvatures, such as those based on normal cycles (see, e.g., 
Cohen-Steiner and Morvan, 2006) or those based on the so-called cotangent formula (see, e.g., Pinkall 
and Polthier, 1993), followed by (c) dividing the resulting integrated curvatures by appropriate area 
terms. The resulting pointwise curvature approximations ki,k2 : V ^ R. are at first only defined on the 
vertex set V of the underlying net. However, we may regard these curvature approximations as functions 
: M — > M. by extending them from V in a piecewise constant manner to the intrinsic Voronoi regions 
of the set V on M. Assuming this extension, we show: 

Theorem 1. Let M be a smooth compact oriented surface without boundary^ immersed into E^. Con- 
sider a discrete net of curvature lines on M such that at each vertex the sampling condition (10) is 
satisfied. Let e be an upper bound for the edge lengths of the net such that additionally the intrinsic 
e-balls around vertices cover all ofM. Then 

sup \ki{p) - Ki{p)\ <Ce, i = 1,2 , 

peM 

where C depends only on properties of M and the shape regularity (9) of the net of curvature lines. 
A few remarks seem pertinent before proceeding: 

• Uniform pointwise convergence of principal curvatures obtained by a 1-local construction from nets 
of curvature lines is somewhat surprising since in general 1-local polyhedral curvatures may not even 
converge in L^, even if the mesh vertices reside on the smooth limit surface (see Hildebrandt et al., 
2006). 

• Our pointwise curvature approximations arise from dividing integrated "Steiner-type" curvatures by 
associated area terms. For example, the integrated mean curvature of an interior edge of a polyhedral 
surface may be defined as the product between the length of that edge and the signed angle between 
the normals of its adjacent faces - a definition that arises from Steiner's view of considering offset 
surfaces. To obtain pointwise curvature approximations from discrete integrated curvatures, we divide 
by so-called circumcentric areas. While this approach is not new, see, e.g., Desbrun et al. (2005), our 
convergence result may be interpreted as a justification of this construction provided that the edges of 
a polyhedral surface well approximate the principal curvature directions of a smooth limit surface. We 
prove uniform lower and upper bounds for edge-based circumcentric areas that may be of interest in 
their own right. 

• For our result to hold we require the explicit knowledge of positions of the vertices of a net of cur- 
vature lines on a smooth surface as well as the combinatorics of this net. More precisely, our curvature 
approximations at a vertex p require the position of p and the positions of its direct neighbors (with 
respect to the combinatorics of the net). Note that we do not require the knowledge of the entire net, 

' Surfaces with nonempty boundary can be treated with minor technical modifications. 
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Figure 2: The three generic patterns of curvature lines near an umbilical point, called lemon, star, and 
monstar by Berry and Hannay (1977). 



though. (Given such an entire smooth net, it would be trivial to compute the principal curvatures at its 
vertices.) It would be desirable to drop from our approach the requirement of the exact knowledge of 
vertex positions of a smooth net of curvature lines. Here, one avenue for further study might be to con- 
sider discrete analogues of curvature line nets: so-called principle contact element nets, see Bobenko 
and Suris (2008). 

Our uniform convergence result given in Theorem 1 is a consequence of a corresponding local 
error estimate given in Theorem 2. This error estimate holds up to and including umbilical points, 
where singularities in the curvature line pattern arise. Using a refinement sequence for each of the three 
surfaces shown in Figure 2, we observed numerically that although the shape regularity of the net may 
blow up near umbilics, linear convergence with respect to the maximum edge length e remains valid in 
these cases. Our experiments also indicate that linear convergence is optimal. 

Alternative approaches An alternative point of departure for establishing pointwise convergence of 
discrete curvatures is to give up ^-locality and to allow for ^ ^ oo as the mesh refinement increases. 
In fact, the above mentioned convergence results of (ii) and (iii) may be interpreted in this way: by 
decreasing the diameter of the domains over which discrete curvatures are integrated (measured), while 
simultaneously increasing the mesh refinement inside these domains at a sufficiently fast rate, one re- 
covers classical pointwise notions of smooth curvatures in the limit. In a similar fashion, Belkin et al. 
(2008) proposed a discrete Laplace operator based on the heat kernel. This operator converges in a 
pointwise manner if the kernel is scaled down while the mesh resolution is increased sufficiently fast 
relative to the scaling of the kernel. In contrast to these works, which need to allow for ^ oo to 
establish pointwise convergence, our result is obtained by working with the simplest and most local 
definition: k = \. 

2 Discrete curvatures from nets of curvature lines 

In order to motivate our definition of discrete curvatures for nets of curvature lines, we recall some 
important notions of curvature for polygonal curves and polyhedral surfaces. For a similar discussion, 
we refer to Sullivan (2008). 

2.1 Discrete curvatures of "Steiner-type" 

Integrated curvatures for polygonal curves Generalizations of classical smooth notions of curvature 
date at least back to Steiner (1840), who considered parallel offsets of convex hypersurfaces, relating 
integrated or total curvatures to changes in length, area, and enclosed volume. For example, for a convex 
curve y c E^, one of Steiner's formulas reads 




y 



(1) 
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Figure 3: Applying Steiner's formula (1) to the three depicted definitions of offset curves for a given 
polygonal curve lead to the three discrete curvatures in (2). 

where I is the length functional, k denotes the curve's curvature, and is the offset curve obtained by 
displacing y along its normals by some constant amount e. 

Steiner's offset formula can be extended to the non-smooth and non-convex case (Federer, 1959; 
Wintgen, 1982; Zahle, 1986). In particular, various notions for curvatures of polygonal curves may be 
interpreted using Steiner's framework. Consider, e.g., 

kpel^Op, 2sin^, 2tan^| , (2) 

where p denotes an inner vertex of a polygonal curve, and 9p is the turning angle between the two line 
segments incident to p. These notions arise by applying (1) to the three different types of offsets depicted 
in Figure 3. Among these, the first notion is the one considered by Steiner, the second corresponds to 
a finite element discretization using piecewise linear functions, and the third also arises in the theory of 
discrete integrable systems (Bobenko and Suris, 1999; Hoffmann, 2008). 

Integrated curvatures for polyhedral surfaces By a polyhedral surface, we mean a piecewise lin- 
ear immersion of a compact simplicial surface into E^*. Extending the notions of discrete curvatures 
from polygonal curves to oriented polyhedral surfaces leads to the following edge-based definitions of 
integrated normal curvature: 

k^ei^eMl 2sin|||e||, 2tan|||e||} . (3) 

Here 6e £ {-^, ^) is the signed angle between the normals of the two flat faces incident to the edge e. 
Notice that measures curvature orthogonal to e, since there is no curvature along e itself. Integrated 
mean curvature is accordingly defined as //e - 

In the planar limit (^e 0), the definitions in (3) agree up to second order in the angle variable. 
Therefore, as it turns out, it suffices for our purposes to prove convergence of one of these definitions 
in order to obtain convergence for all of them. Convergence of the first definition in (3) in the sense of 
measures was investigated in Fu (1993); Cohen-Steiner and Morvan (2006). 

For completeness, we remark that the above edge -based definitions give rise to vertex-based notions 
of integrated mean curvatures by adding the mean curvatures over all edges emanating from a given 
vertex p, i.e., 

Hp=^-J]H,. (4) 

The factor i takes the meaning of distributing the normal curvature of each edge equally among its two 
adjacent vertices. 

Finally, the scalar-valued definitions considered so far can be extended to corresponding vector- 
valued notions. In the edge-based case, we obtain normal curvature vectors kg by multiplying with 
the angle-bisecting unit normal vector at e (see Figure 4, left), and similarly for mean curvatures. Anal- 
ogously to (4), we then obtain vertex-based mean curvature vectors. We remark that for - 2 sin y ||e||, 
the resulting mean curvature vector coincides with the surface area gradient at p when restricting to 
piecewise linear surface variations (yielding the so-called cotangent formula, see Pinkall and Polthier, 
1993). Its convergence in the sense of Sobolev norms was studied in Hildebrandt et al. (2006). 



4 



Figure 4: Edge-based quantities. Left: dihedral angle 9^ and discrete curvature vector kg. Right: dual 
edge and circumcentric area Ag. 



From integrated to pointwise curvatures In order to obtain pointwise curvatures, we divide the 
above integrated curvatures by corresponding area terms. Intuitively, these areas can be thought of as 
the domain of integration from which integrated curvatures were obtained. Whether or not one obtains 
convergent curvatures depends on a careful choice of these areas. It turns out that for triangulated poly- 
hedral surfaces, one good choice are the so-called circumcentric areas, such as considered in Desbrun 
et al. (2005). For each edge e we define 



where ||*e|| denotes the intrinsic length of the circumcentric dual edge ★e. This dual edge intrinsically 
connects the circumcenters C\ and C2 of the two triangles T\ and T2 incident to e. Here, intrinsic means 
that one can think of T\ and T2 as being unfolded onto the plane (see Figure 4, right). The sign sgn(*e) 
is positive if along the direction of the ray from Ci through C2, triangle T\ lies before T2, and negative 
otherwise. Note that sgn(*e) < (and therefore Ae < 0) iff ae + y^e > n, where and y6e are the 
angles opposite to e in the triangulation (see Figure 4, right). Consequently, we require lower bounds 
that ensure positivity of circumcentric areas. For nets of curvature lines, we provide such bounds in 
Section 3.1. 

Similar to vertex-based integrated curvatures, we obtain vertex-based circumcentric areas from the 
edge-based case via 



where the sum is taken over all edges emanating from a given vertex p. If all edges e incident to a 
vertex p are intrinsically Delaunay (compare Bobenko and Springborn (2007)), then Ap coincides with 
the intrinsic Voronoi area of p and is therefore positive. However, as pointed out in (Dyer and Schaefer, 
2009), Ap might become negative in general, and a bound similar to the edge -based case is not possible. 
Therefore, we will not treat vertex-based pointwise curvatures based on Ap. 

2.2 A local error estimate for discrete curvatures 

In this section, we state our main local error estimate (Theorem 2), from which we derive our global 
uniform convergence result (Theorem 1). 

Throughout we assume that M is a smooth compact oriented surface without boundary immersed 
into E^. By a discrete net on M we mean a cellular decomposition of M such that all attaching maps 
are homeomorphisms and the intersection of any two cells is either empty or a single cell. As usual, we 
denote by E the set of edges and by V the set of vertices. We also assume that all edges are smoothly 
embedded. In a discrete net of curvature lines on M, all edges are additionally required to be segments 
of curvature lines, non-umbilical vertices are required to have valence four, and umbilical vertices are 




(5) 




(6) 



5 



Figure 5: The triangulated vertex star at a vertex /? is a polyhedral surface approximating the (curved) 
discrete net of curvature lines at p. 



required to have valence greater than two. In a completely umbilical region (such as S^), any net in the 
above sense serves as a net of curvature lines for our purposes. 

In order to be able to apply the concepts of discrete curvatures on polyhedral surfaces to nets of 
curvature lines, we require local polyhedral approximations of smooth curvature line nets. 

Local polyhedral approximation In the sequel, the letters e, f, and g will be reserved for (curved) 
edges in the edge set E, incident to a common vertex p, while corresponding bold face letters, e, g, 
and f , will denote the straight edge vectors in obtained by connecting the endpoints of e, f, and 
g, respectively, by straight lines (see Figure 5). We additionally assume that these edge vectors are 
oriented such that they point away from p. For each disjoint edge pair {e, f) incident to p and contained 
in a common 2-cell, we consider the flat triangle spanned by e and f. The union of these triangles 
forms the triangulated vertex star of p, denoted by st{p). Whenever we consider the triple (e, f , g), we 
will always assume that the pairs (e, f ) and (e, g) span two triangles in st{p), such that e becomes their 
common edge. Finally, as later justified by our sampling condition (10) and Corollary 9, we may assume 
that n • (e X f ) > and n • (e x g) < 0, where n denotes the normal of M at p. 

Each triangulated vertex star st{p) thus yields the requisite local polyhedral approximation^, which 
forms the basis for our curvature approximations. As outlined in the previous section, our definition of 
pointwise curvatures relies on the division by certain circumcentric areas, which may become zero or 
negative in general. This motivates, for a given principal direction, to choose the associated edge vector 
with maximal circumcentric area. 

Definition 1 (area maximizing edge). Consider a vertex p in a discrete net of curvature lines. If p 
is umbilical, we call an edge vector e area maximizing ;/ it maximizes the circumcentric area among 
all edges emanating from p in the local polyhedral approximation. If p is non-umbilical, let vi be 
the principal direction canonically associated with an edge vector e. We call e area maximizing if it 
maximizes the circumcentric area among the two edge vectors associated with vi. 

We show in Section 3. 1 that area maximizing edges always have circumcentric areas that are bounded 
away from zero. 

Definition 2 (principal curvature approximations). Consider a vertex p in a discrete net of curvature 
lines and let e be an area maximizing edge (associated with a principal direction vi ;/ p is non- 
umbilical). Then 



klip) ■■= 




^Observe that we do not require that our local polyhedral approximations yield a consistent global one. 
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Figure 6: A discrete net of curvature lines, with the two families of directions depicted as solid and 
dotted lines, respectively, together with the circumcentric areas (gray and white diamonds) per edge. 



defines the principal curvature approximation of K2{p), where K2 refers to the principal curvature corre- 
sponding to V2 if p is non-umbilical and refers to the unique normal curvature if p is umbilical. Here kg 
is one of the edge-based integrated polyhedral curvatures defined in (3). 

The fact that in the above definition the edge vector e is associated with vi while ^ approximates 
K2 is not an oversight: measures curvature orthogonal to e. 

The intuitive reason for dividing by twice the circumcentric area in the above definition may (at 
least qualitatively) be explained as follows. Consider Figure 6, where one set of principal curvature 
directions, say those associated with vi, is depicted by solid lines, while the direction corresponding to 
\2 is represented by dotted ones. Likewise, the circumcentric areas corresponding to vi -directions are 
drawn as gray diamonds, while the circumcentric areas corresponding to V2-directions are represented 
as withe diamonds. Roughly, the gray diamonds cover only half of the total surface. Hence, taking 
only the gray diamonds as regions of support for our (integrated) principal curvature approximations 
corresponding to ki, would mean to be roughly missing a factor of two. This motivates, for each edge 
along a given principal direction, to consider twice its circumcentric area as the domain of integration. 

Global constants For each p e M, let S (p) denote the shape (or Weingarten) operator. Our estimates 
depend on both S and its covariant derivative, . Accordingly, we define 

'K ■- max \\S{p)\\op and 7C' ■- max (max ||Vv5(p)||op | , (7) 

P P \llv||=l / 

where || • ||op denotes the usual norm for linear operators. Note that 7C is an upper bound for the normal 
curvatures of M, whereas "TC' provides an upper bound for directional derivatives Vy(/c,) of principal 
curvatures /c,-. 

Local constants We also consider local constants - shape regularity p and maximum edge length e - 
that are specific for each vertex in the net of curvature lines. The reason for introducing local constants 
is that a high aspect ratio at one vertex should not affect the sampling condition (see below) at another 
vertex. In the following, we assume an arbitrary but fixed vertex p. 

We let e denote the largest intrinsic edge length over all (curved) edges e € E emanating from p eV 
(denoted by e ~ p), 

€ = max 1(e) . (8) 

Notice that the length of every edge vector emanating from p in st{p) is thus also bounded above by e. 

Our estimates also depend on shape regularity, or aspect ratio. We define p > 1 to be the smallest 
number such that for all pairs (e, f ) of edge vectors emanating from p and forming a triangle in st{p), 
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one has 



- < llell as well as < p . (9) 

P lie X f|| 

The former inequality implies ^ < ||||| < p, while the latter means sin Z(e, f) > ^. 

Sampling condition In addition to the above definitions of maximum (local) edge length and (local) 
shape regularity, we assume the (local) sampling condition 

In some of our estimates, it will suffice to work with weaker sampling conditions, such as e < or 
^ - TK' ^^^^ of which are implied by (10). 

Theorem 2 (Local error estimate). Let M be a smooth compact oriented surface without boundary 
immersed into E^. Consider a vertex p € V in a discrete net of curvature lines on M and assume the 
sampling condition (10). Let {ki{p), k2{p)) denote the approximations of the smooth principal curvature 
{K\{p),K2{p)) as in Definition 2. Then 

Mp)-Ki{p)\<Ce i=\,2. (11) 

The constant C = C^K, "K' ,p) depends only on the curvature bounds (7) and the shape regularity (9). 

Remark. Equivalent estimates can be obtained when replacing the scalar- valued quantities in (11) by 
their corresponding vector-valued counterparts (for definitions, see Section 2). Proofs remain nearly 
identical. 

Our global uniform convergence theorem stated in the introduction is a direct consequence of our 
local error estimate. 

Proof of Theorem 1. Let / : M — > V be defined by mapping each p e M to its nearest point in the 
vertex set V, i.e., dM{p, f{p)) - duip, V), where dM denotes the intrinsic distance on M. We extend 
our curvature approximations (initially only defined on V) to functions fc, : M ^ R via p i-^ ki{f{p)). 
By the assumptions of Theorem 1, e was globally chosen such that dMip, fip)) < e for all p € M. This 
implies, by connecting p and f(p) by a shortest geodesic arc, that 

\Kiip) - Ki{f{p))\ < nc'e . 

Furthermore, e was globally chosen such that it maximizes global edge length. Hence Theorem 2 
implies that there exists a constant C such that 

\h{f{p)) - Ki{f{p))\ < Ce 

for all p e M. Additionally, note that the constant C in Theorem 2, besides depending on "TC and "TC', 
is monotonically increasing with respect to shape regularity. (This will become evident in the proof of 
Theorem 2.) Hence, C depends only on 7C, 7C', and the largest local shape regularity constant p. This 
together with an application of the triangle inequality, 

Mp) - Ki{p)\ = \kiifip)) - Ki{p)\ < Mfip)) - Ki{f{p))\ + \Kiifip)) - Ki{p)\ , 

implies the claim. □ 
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3 Proof of local error estimate 



The proof of Theorem 2 proceeds in several steps. First, we provide uniform lower and upper bounds 
for the edge-based circumcentric areas by which we divide integrated curvatures to obtain pointwise 
notions (Section 3.1). In a second step, we provide estimates for edge-based integrated curvatures by 
using their corresponding discrete curvature vector. We establish that for each vertex p in a. net of 
curvature lines, the projection of these vectors onto the tangent plane TpM is negligible. Furthermore, 
we show that the remaining normal component leads to the error estimate in Theorem 2 up to a certain 
error term (Section 3.2). While for general meshes the resulting error term cannot be controlled (and 
indeed causes failure of pointwise convergence), we provide bounds for this error term for the specific 
case of nets of curvature lines (Section 3.3). 

Basic assumptions In order to avoid excessive repetition, we summarize our basic assumptions and 
notations. We write p for a non-boundary vertex of a polyhedral surface, with vertex star denoted by 
st{p). We assume that e is an upper bound for the length of the (straight) edges emanating from p, 
and we let p denote the shape regularity as defined in (9). If st{p) arises from the local polyhedral 
approximation of a net of curvature lines, then e is defined by (8), i.e., as the maximum edge length of 
the curved edges emanating from p. Throughout, we assume the sampling condition (10). As before, 
we denote curved edges by e, f, and g, and their corresponding straight edge vectors by e, f , and g. 

3.1 Uniform bounds for circumcentric areas 

In this section we prove upper and lower bounds for the circumcentric areas of area maximizing edges 
in the sense of Definition 1 . 

Proposition 3. Consider a vertex p in a discrete net of curvature lines and let e be an area maximizing 
edge vector at p. Then our basic assumptions imply the existence of some C > such that 

where C only depends on the shape regularity constant p. 

Note that for non-umbilical vertices, this result implies the existence of an edge with positive cir- 
cumcentric area for each of the two principal curvature directions. 

The remainder of this section is concerned with proving Proposition 3. First observe that 

Ae = ^(cotae + cot^e)lle|P - /"^''^ Help , 
4 4 sm Qe sm pe 

where ae and ySe are the angles opposing e in the two triangles meeting at e, respectively. The requisite 
upper bound on Ae is relatively straightforward to obtain. 

Lemma 4 (upper bound). Let p be a vertex in a discrete net of curvature lines on M. Then our basic 
assumptions imply that each edge e emanating from p satisfies 



Proof. Clearly, we have 



A.<pV 



A. . 



4 sin ae sin fie 

Let f be the edge vector emanating from p such that ae belongs to the triangle formed by e and f . Then 
the definition of shape regularity (9) implies 

. ^ ||(f-e)xf|| ^ llexfll ^ 1 ^ J_ 
||(f - e)|| ||f II - (p + l)||e|| ||f II - (p + l)p - 2p2 • 

A similar estimate holds for Pe- Hence, Ae < p'^llelP < p^e^. □ 
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Similarly, we obtain a lower bound for at least one edge emanating from p. 



Lemma 5 (lower bound). Let p be a vertex in a discrete net of curvature lines on M. Then our basic 
assumptions imply that there exists an edge vector e emanating from p such that 



1 

4p3 



,2 



Proof. Let e be the shortest edge emanating from p, let f be the straight edge emanating from p such 
that Q-g belongs to the triangle formed by e and f, and define y :- Z(e, f). Since ||e|| < ||f ||, it follows that 
Iffe ^ - y)- Since in particular < ag < f , it follows that cotae > f - By (9) we have siny > ^, 
and hence 

n y siny 1 

COtc.e>--ae>-> — >-. 

Applying similar arguments, we obtain cot/Je > j^- Together, this yields 

Ae - ^(cottte + cot;Se)lle||^ > ^M? > -^e^ . □ 
4 4p 4p^ 



Lower bounds for vertices of valence four The above lower bound on the circumcentric area of at 
least one edge suffices for umbilical vertices. However, it does not suffice at non-umbilical ones, since 
we require lower bounds for edges associated with each of the two principal directions. We achieve this 
by showing that for each vertex of valence four, there are at least three edges that satisfy the required 
lower bound. Hence, for each vertex in a net of curvature lines, we have at least one good edge per 
principal direction. 

We note that some of the following results are also valid for vertices of valence different from four, 
such as, in particular. Corollary 9. 

As before, we let and j6e denote the angles opposing the straight edge e in the two triangles 
meeting at e, respectively. If ae > 5, /3e ^ 3, and + fie < n - 6 for some 6 > 0, then 

A, . , (13) 

which provides a useful lower bound if 6 can be bounded away from zero. Accordingly, we introduce 
the notion of 5-Delaunay edges, a nomenclature that is borrowed from the classical case of Delaunay 
triangulations (corresponding to 5 = 0). 

Definition 3 (5-Delaunay). Let Oe and /3e be the angles opposing an edge e in the two triangles meeting 
at e, respectively. Then e is called 6-Delaunay if there exists 5 > Q such that Oe > 6, fie ^ ^, cind 
ae+ Pe<^ - 5. 

Assume for a moment that p has valence four and that st{p) is planar. Assume further that all of 
the eight angles opposing the four edges emanating from p are bounded from below by 5. Then it is 
straightforward to verify that at least three among the four edges incident to p are 5-Delaunay. In general, 
however, st{p) is not planar, and we have to account for Gaussian curvature. As usual, we define discrete 
Gauss curvature at a vertex p ofn polyhedral surface as the angle defect, i.e., by Kp = 2n- y,-, where 
y, are the intrinsic angles meeting at p. We obtain: 

Lemma 6. Let p be a non-boundary vertex of valence four on a triangulated polyhedral surface, and 
let ai,(ij denote the pairs of angles opposing the four edges emanating from p. Assume that a,- > S and 
/?( > Sfor all i - I,. . . ,4, as well as Kp < 25. Then at least three among the four edges meeting at p 
are 6-Delaunay. 

Proof. Using Kp = -In + ZiLiCo^i + Pd^ the result follows from a straightforward calculation. □ 
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We now show that our basic assumptions imply the assumptions of Lemma 6 when setting S '■- 
To see this, first observe that with the above notations, ae > sinae > 6hy (12), and analogously for y&e- 
It remains to check the condition Kp < 26 for the discrete Gauss curvature. The requisite bound will 
be established in Lemmas 7 and 8. The resulting consequence for the existence of 5-Delaunay edges is 
summarized in Lemma 10. Finally, Lemma 11 establishes the lower bound for A^. 

Lemma 7. Let p be a non-boundary vertex of an oriented polyhedral surface. Assume that the normals 
of the triangles incident to p make an angle no greater than (p € [0, tt) with some fixed direction in B?. 
Then the discrete Gauss curvature associated with p satisfies Kp < 2n{\ — coscp). 

Proof. The lemma is a consequence of the Gauss-Bonnet theorem. Let iii,n2, . . . ,n/ denote the unit 
normals of the triangles T\,T2, . . . ,Tj incident to p, ordered according to the orientation of the polyhe- 
dral surface. Each n, represents a point on S^. Connecting consecutive pairs (n,, n,+i) by geodesic arcs 
on yields a spherical polygon P (possibly with intersecting edges). For each /, the exterior angle of 
P at n, is then equal to the interior angle of the Euclidean triangle T, at p. In particular, the sum E of the 
exterior angles of P satisfies 1. - 2n - Kp. 

Moreover, by assumption we have € [0, tt), so all n, lie on the same hemisphere. We can hence 
consider the spherical convex hull P of ni, n2, . . . , n;, i.e., the smallest spherical polygon that contains 
all n, and that is convex with respect to shortest geodesic arcs. Let S denote the sum of exterior angles 
of P. It is easy to verify that E < E. Hence Kp <2n -f. - area{P), where the last equality follows from 
the Gauss-Bonnet theorem. Since (p € [0, n), the polygon P is contained in a geodesic disk of radius (p, 
the area of which is 2n{\ - cos 0). This proves the claim. □ 

In order to make use of the previous lemma, we seek a bound on the angle cp between the surface 
normal at p e M and the normals to the triangles incident to p. Note that our basic assumptions, and in 
particular our sampling condition, imply f < where "7C is our curvature bound. Hence we can infer 
from Morvan and Thibert (2004, Section 3, Corollary 1): 

Lemma 8. Let p be a vertex (umbilical or not) in a net of curvature lines on M. Given our basic 
assumptions, st{p) can be oriented such that the maximum angle cp € [0, n) between the surface normal 
at p e M and the normals to the triangles in st{p) satisfies sin (p < (4p -I- 2) TCe, where p is the shape 
regularity. In particular, sin < |- 

Corollary 9. Under the assumptions of Lemma 8, st{p) can be oriented such that the orthogonal pro- 
jection of st{p) onto TpM is infective and orientation-preserving. 

Note that our assumptions are slightly different from those used in Morvan and Thibert (2004); in 
fact, our assumptions are stricter. While their sampling condition bounds the maximal extrinsic distance 
between two neighboring vertices, we consider the intrinsic length on the smooth surface, which is 
always larger. Additionally, Morvan and Thibert (2004) require that the distance between the discrete 
and the smooth surface is less than the reach of the smooth one. This requirement is implicitly fulfilled 
locally by the sampling condition e < since the reach of a surface patch formed by an intrinsic 
6-disk around p is nothing but the minimal radius of curvature of that surface patch . 

Lemma 10. In addition to the assumptions of Lemma 8, assume that p is of valence four Let 5 :- 
where p is the shape regularity. Then at least three edges incident to p are 6-Delaunay. 

Proof. Using that (1 - cos0) < sin^0 for (p e [0, |], Lemmas 7 and 8 show that the discrete Gauss 
curvature at p satisfies 

Kp < 2n{Ap + 2f'K^e^ < {\6'Kpef . 

Setting 6 = j^, the sampling condition (10) implies Kp < 25. Moreover, (12) implies aj > 6 and 
Pi > 5 for all pairs of angles opposing the four straight edges emanating from p. Finally, Lemma 6 
implies the claim. □ 
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Lemma 11 (lower bound for valence four). Under the assumptions of Lemma 10, there exist at least 
three straight edges among the four edges emanating from p, such that 

16p4 

for each edge e among these three. 

Proof. Observe that (13) implies 4Ae > sin5||e|p for all (5-Delaunay edges. Applying Lemma 10 and 
using sin > ^ gives 



Ae > ; sini ^ I ||e||2 > ^ > -^e^ 



16p2 - I6p4 



Proof of Proposition 3. Lemma 4 provides an upper bound for both umbilical and non-umbilical ver- 
tices. Lemma 5 provides the requisite lower bound for umbilical vertices. Finally, Lemma 1 1 provides 
the lower bound for non-umbilical ones, since it implies the existence of at least one edge per principal 
direction with circumcentric area bounded from below. □ 



3.2 Estimates for discrete integrated curvatures 

In this section, we establish a bound for the difference between the edge-based integrated curvatures ke 
and (an appropriately scaled version of) the smooth principal curvatures of M. Specifically, for a given 
edge e in our polyhedral approximation, we work with the discrete integrated curvature kg - 2 sin Y||e|| 
introduced in Section 2. 1 . 



Darboux frames For our purposes, it turns out to be useful to express vectors in frames that are locally 
adapted to the geometry of the surface M. Specifically, a Darboux frame at a non-umbilical point p e M 
is an adapted frame given by (vi,V2,n), where vi and V2 are (normalized) principal directions of M at 
p, and n = vi X V2 is the surface normal (induced by the orientation of M). For umbilical points, any 
adapted (i.e., Vi,V2 e TpM) and orthonormal (i.e., ||vi|| = ||V2|| = 1 and n = vi x V2) frame may be 
considered a Darboux frame. Throughout, we employ the notation 

e = (ei,e2,e„) 

to represent a vector e in the coordinates given by a Darboux frame. In the sequel, we assume a fixed 
Darboux frame at every vertex p of our discrete net of curvature lines. If p is non-umbilical then each 
edge vector emanating from p is canonically associated to exactly one of the principal directions vi 
or V2. If p is umbilical, we additionally require an explicit association of each edge vector with one 
of either vi or V2. In order to state the main result of this section, we require the notion of tangential 
deviation of an edge vector with respect to a Darboux frame. 

Definition 4 (tangential deviation). Let p be a vertex in a discrete net of curvature lines on M and let 
(vi , V2, n) be the Darboux frame at p. If an edge vector e is associated with the principal direction y\, we 
call Cd - \e2\ its tangential deviation, see Figure 7. Likewise, ife is associated with \2, we let = \ei\. 

The following proposition summarizes the main result of this section. 

Proposition 12. Using the notations of Definition 4, let the edge vector e emanating from vertex p be 
associated with vi, and let its directly neighboring edges f and g be associated with V2. Then our basic 
assumptions imply that 

\ke - 2A,K2\ < C [\6,(ed + fd + gd)\ e + e^) , 

where K2 denotes the principal curvature ofM in direction V2, 6^ = K2 — k\ (with 6k = in the umbilical 
case), and C = C(7C, 7<"',p). 
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'p + Rvi 



Figure 7: The tangential deviation e^ of a straight edge e. 



Remark. We point out that the estimates in this section are not entirely specific for nets of curvature lines. 
In fact, they hold in the more general setting of arbitrary smooth nets embedded in M, as long as we 
require bounded shape regularity and the sampling condition e < where e is the maximum intrinsic 
edge length of the net and 7C denotes our usual curvature bound. As a consequence, the estimates 
of the current section do not suffice to establish our main error estimate. Indeed, to obtain uniform 
convergence, we require that the term d^ied+fd+gd) appearing in Proposition 12 is of order e^. However, 
this is false for general nets (and causes failure of uniform convergence), but is true for nets of curvature 
lines as we will show in Section 3.3. For clarity's sake, though, we have decided to restrict the discussion 
of the current section to nets of curvature Unes and rely on our (rather strong) basic assumptions set forth 
in the beginning of Section 3. 

Proposition 12 is proven in several steps. We commence by estimating the normal component e„ of 
the straight edge vector e (Lemma 13 and Corollary 14). Then we switch to the edge-based curvature 
vector kg corresponding to kg. We show that its tangential component is negligible (Lemma 15). The fi- 
nal proof of Proposition 12 is given at the end of this section, where we show that the normal component 
of ke yields the desired estimate. 

Lemma 13. Let p be a vertex of a discrete net of curvature lines on M. Consider an edge e e E 
emanating from p with corresponding straight edge vector e. Writing e - {e\,e2,en) with respect to a 
Darboux frame centered at p, our basic assumptions imply that 

\en\ < TCe^ < TCe^ . 

Furthermore, if P{x,y) = '^x^ + denotes the osculating paraboloid, then 

\en-P{eue2)\<Ce\ 

where C only depends on our global curvature bound and the bound on curvature derivatives 7C'. 

Proof. Using a Darboux frame at p, the surface M can locally be parameterized by a height function 
h{x,y) over the tangent plane TpM. In the coordinates of the Darboux frame, we have e = {e\,e2,en), 
with en = h{e\,e2). We let d = ||(ei,e2)ll> where throughout this proof || • || denotes the Euclidean norm 
in the parameter domain. Furthermore, we consider the constant vector field y{x,y) - {ei,e2)/d in the 
parameter domain. 

Let D'^h denote the ith iteration of the directional derivative of h along v, i.e., D^/j = Dyh and 
- Dy{D'~^h), and let D'h denote the ith total derivative with respect to the standard Euclidean metric 
in the parameter domain. Observe that h{0) = Dh{0) = 0, and hence 

en^h{eue2)^ f f Dlh{T\)dTdt . (14) 
Jo Jo 
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Consequently, in order to prove the first part of the lemma, we seek an upper bound on in terms of 
7C. 

Let S denote the shape operator and let !(•, •) and ][(•, •) = 1{S •, •) be the first and second fundamental 
form of M, respectively, with respect to the local parameterization induced by h. From 



and 



we obtain the estimate 



I(u, v) = u^(Id + Dh Dh^)y (15) 



I(v,v) 1(5 v,v) D^h 

Ky = = = , (16) 

I(v,v) I(v,v) Vl + ||D/i||2l(v,v) 



\Dlh\ = \Ky\ Vl + mmy, V) < 7C(1 + ||D/j||2)i . (17) 



In order to bound \\Dh\\, first observe that our sampling condition implies e < This, in turn, provides 
a bound on the (positive) angle between the surface normal iip at p and the surface normal at any 
point q on the (curved) edge e (Z M incident to p in the given net of curvature lines on M. To see this, 
consider an arc-length parametrized curve y : [0,^] ^ M with 7(0) = p and y(^) = q. Notice that our 
basic assumptions imply that y can be chosen such that ^ < e. The Gauss image -y = noyofyisa curve 
on the unit sphere. The length of y is therefore bounded from below by /(n^, n^), i.e., the length of the 
minimizing geodesic joining and on S^. The tangent vector of y at a point y{s), s e [0, ^], is given 
by Sy'{s), and the norm of this vector is therefore bounded above by 7C. Hence, 



1 

4np,ii,)< \\Syis)\\ds<^^<^e<-. 



Writing q - (qi,q2, qn) with respect to the Darboux frame centered at p - (0, 0, 0), it follows that 

l|£>%i,^2)ll - tanz(np,iig) < tan ^ . (18) 

Plugging this into (17), we find that |D^/j| < I'K, which, together with (14), yields the first part of the 
lemma. 

In order to prove the second part, we first note that D^h{0) = D^P{0) and that D^P is constant. 
Hence 



IKeuei) - P(ei,e2)\ = 



r r fo 

Jo Jo Jo 



^h{o-y) dcrdTdt 



(19) 



We consequently seek a bound for \Dyh\ in terms of 7C and 7C'. Considering (16) and taking another 
derivative with respect to v yields 

Dlh = Dy(Ky) Vl + \\Dhfl{\, V) 

+ KyDy[yl\ + \\Dhf)l{y,y) (20) 
+ KyyJ\ + \\DhfDy{l{y,y)). 

We bound the terms appearing on the right hand side one by one. To treat the first term, we use (16) and 
V/ = 0, where V denotes covariant differentiation with respect to the metric induced by I, to derive 

I((Vv5)v, V) + 21(5 V, Vyv) - iKyljy, V,v) 

DyKy ^ . (21) 

I(v,v) 
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Figure 8: The discrete curvature vector ke is the sum of JfC (the rotation of e by | about the axis e X f ) 
and JgC (the rotation of e by | about the axis e x g). 



Let u(x, 3^) = u be a constant vector field in the parameter domain. Using the Koszul formula for the 
Levi-Civita connection and applying (15) yields 

2I(VvV,u) - 2v(I(v,u)) - u(I(v,v)) = IDlhD^h . 

The last equality only depends on the value of u(x, y) at the point (x, y) and therefore holds for any field 
u in the parameter domain. From (17), (18), and ||v|| = 1, we obtain 

|I(VvV,u)|<C||u||, 

with C = CCTC). This can be used in (21), together with our bounds on the norms of S and VS in terms 
of "TC and respectively, to obtain an upper bound for the first term in (20) by C(7C,'7C')- Bounds 
for the remaining two terms in (20) can be obtained in a similar fashion using (15) and (16), proving 
the estimate \Dlh\ < C(7<",7("')- Using (19) and d < e then implies the claim of the second part in the 
statement of the lemma. □ 

Corollary 14. With the same assumptions as in Lemma 13 and by defining = K2- k\, we obtain 

.„^|e2 + |.^0(.^)^|e2-|e? + 0(.^), 
with\0{e')\ < C{'K,'K')e'. 

Proof. From the second part of Lemma 13 and simple algebraic manipulations, we deduce that 

e„ = - '^el + %el + Oie') = ^e^ + ^^.^ + ^(^3^ _ 
2 2 2^ 2 2^ 

The last equality follows from Lemma 13 and the sampling condition e < (which is implied by our 
basic assumptions), which together ensure that e^ < TC^e^ < ye^. The second equation in the statement 
of the corollary follows analogously. □ 

For the following discussion, it will be useful to work with the curvature vector ke corresponding to 
the discrete integrated curvature kg = 2 sin Yl|e|| (see Section 2.1). With our usual notions, consider the 
two edge vectors f and g directly neighboring e in st{p) (see Figure 8). Recall from Corollary 9 that we 
may assume that (e x f ) • n > and (e x g) • n < 0, where n is the normal of M at p. A straightforward 
calculation reveals that we can express the discrete curvature vector by 

ke ^ JfC + JgC , (22) 
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where Jf and Jg denote rotations by | around the axes e x f and e x g, respectively (see Figure 8). 
Consider now the splitting 

ke = (ke); + (ke)„ 

of ke into its tangential and normal component with respect to TpM. We show that the tangential 
component is negligible: 

Lemma 15. Under the same assumptions as in Lemma 13, the projection of the discrete curvature 
vector ke onto the tangent plane TpM satisfies ||(ke)r|| < lOTC^p^ e^. 

Proof. With the notions and assumption of the preceding discussion, we first note that a straightforward 
calculation reveals that 

e"f-(f-e)e 

and analogously for JgC. 

With respect to TpM, let Cf, f,, and gt denote the tangential components of the vectors e, f, and g, 
respectively. Then, by assumption, we have (e? x fj) • n > and (e? x gj) • n < 0. The tangential part of 
the discrete curvature vector ke is given by (ke)r - (JgC), + (Jff)?, where we deduce from (23), using the 
coordinates of our fixed Darboux frame at p, that 

||e|pf,-(f-e)e, ||e,xf,|| ^ el^-enfn^t 
(Jfe> = = ■^{-e2,e]_,Q) + 



llexfll llexfll llexfll 

and similarly 

,T X llerXg^ll elgt-eng„tt 

(Jge)f = r-(-e2,ei,0) + — — . 

Ilexgll llexgll 

Using the first part of Lemma 13 and the definition of shape regularity (9), we obtain 

lie X f||2 - lie, X f,||2 + {ej, - ejnf + {ejj - ejfnf 
< lie, X l,t + 8(7C||e|| ||f|k)2 

'2„2^2|| 



Therefore, we have 



< lie, X f,|r + BTCV^IIe x f|| 



llexfll - llexflp 

and similarly for g. It follows that the two terms in (Jfe + Jge), containing {-e2,ei,0) cancel up to a 
term bounded by (87C^p^)6^, where the power of three is due to the fact that the norm of (-^2, ^i, 0) is 
also bounded by e. Moreover, we observe that the first part of Lemma 13 and the definition of shape 
regularity (9) yield 

Welft - enfn^tW ^ "TC^IIeinifll + -TC^HelPHflp 3 



llexfll - llexfll 
and analogously for g. Therefore, we arrive at 

||(ke),|| - IKJfe + Jge),|| < (lOTCV)^' , 

proving the claim. 

The above implies the main result of this section: 
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Proof of Proposition 12. To estimate the normal component (ke)„ - (Jfe)„ + (Jge)„ of the discrete cur- 
vature vector, we use (23) to obtain 

/T ^ e2/«-(f-eK 
^J^^^" = llexfll ' 

and analogously for (JgC),,. We note that the circumcentric area of e can be expressed in a similar manner 
as a sum Ae - A^j + Ae,g, where 

- 4||e X fll ' ' 
and analogously for Ae,g. Applying Corollary 14, we obtain 

(J^^^" ^ ^Ml ^^^^^ 

= 2AefK2 - 6J} 5 J • e ^ — + Oie') , 

■'^ 2||exf|| 2||exf|| ^ ' 

and analogously for (Jge)„, with < C{^,'K' ,p)e^ . Applying the shape regularity condition (9) 

yields 



|(Jfe)„ - 2Ae,f^2| < C i\6,{fl + f • e)| + 6^) , 



and similarly for the difference between (Jge)„ and 2Ae,g/C2. According to the statement of Proposi- 
tion 12, we assume that e is associated with the direction vi, whereas f and g are associated with V2. 
Using our notion of tangential deviation from Definition 4 and noting that |f • e| < \e2 + /ik + 0{e^), we 
arrive at 

|(ke)„ - 2Ae^2l < C{\6,{ed + f, + g,i)\ e + e^) , 
which implies the claim, since (ke)? is of order O(e^) by Lemma 15. □ 



3.3 Estimates for non-umbilical vertices 

As mentioned before, the results of the preceding section are not entirely specific for nets of curvature 
lines but hold for a larger class of smoothly embedded nets. As such, these results do not suffice to 
prove our main error estimate in Theorem 2, due to the failure of uniform convergence of curvature 
approximations constructed in a 1 -local manner for general nets. Indeed, assuming that e is associated 
with the principal direction vi, we seek a bound of the form 

|/te - 2Aei^2l < Ce^ , (24) 

from which we may derive the desired main error estimate by employing our results on the existence of 
uniform lower and upper bounds of (sufficiently many) circumcentric areas (see Section 3.1). 

While a bound of the form (24) does not hold for general smooth nets, it is indeed valid for nets of 
curvature lines. This is a consequence of the results of the preceding section and the fact that for nets of 
curvature lines we have 

\S,{p)ed\ < Ce^ , (25) 

which is trivially satisfied for umbilical vertices and true for non-umbilical ones provided that we asso- 
ciate e with its canonical principal direction. This is precisely the main result of this section. 

Observe that there is a simple case, where (25) is obviously fulfilled: let M be a paraboloid, and let 
p be its apex. Consider the four edge vectors emanating from in a local polyhedral approximation of 
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a net of curvature lines containing /? as a vertex. Then the tangential deviation of each of these edges 
vanishes, so (25) is clearly satisfied. 

The main difference between nets of curvature lines on arbitrary smooth surfaces and the specific 
case of a paraboloid is the fact that curvature lines usually have non-zero geodesic curvature k^. While 
the tangential deviation Cd can always be bounded by Ce^, with C = C{'K,k^), this does not suffice 
for a uniform error bound, since may blow up at umbilical points. Perhaps surprisingly, though, the 
product \6Kip)ed\ can be bounded for nets of curvature lines: 

Proposition 16. Let e be an edge vector emanating from a non-umbilical vertex p in a local polyhedral 
approximation of a discrete net of curvature lines on M. Ift is associated with its canonical principal 
direction, then our basic assumptions imply 

Mp)ed\<{^^+4<K')e\ 

where Cd is the tangential deviation ofe and 'K, 'K' denote our usual bounds on normal curvatures and 
their derivatives, respectively. 

The intuition behind this statement is as follows. Roughly, Cd is proportional to the geodesic curva- 
ture of the curvature line corresponding to e. This geodesic curvature, in turn, is inversely proportional 
to 6k, as the next lemma shows. Therefore, the product d^ed can be uniformly bounded. 

Lemma 17. At any non-umbilical point of M, the geodesic curvature of the principal curvature line 
along vi satisfies 

^ ^" ^ and thus \ < 



where /ci and denote the principal curvatures corresponding to the principal directions Vi and V2, 
respectively. 

Proof. Since vi and V2 are orthonormal eigenvectors of the shape operator S , we have 5 Vi • V2 = 0. We 
use the Frenet formulas Vy,vi - a^V2, VviV2 = -a^vi, VyjVi - /(^V2, and Vv2V2 = -/<|^i well as the 
Codazzi-Mainardi equation ( Vu5 )v - ( VyS' )u to obtain 

0-Vv,(5Vi-V2)- Vv,(5V2-Vi) 

- (Vv,5)v2 • vi + 5(Vv,V2) • vi + 5V2 • Vv.vi 
= (V,,5)vi • vi + 5(-4vi) • vi + 5V2 • (4v2) 



= (Vy,S)Vi • Vi -I- K^A-Kiyi ■ Vi + /f2V2 • V2) 



= Vy^iS Vi • Vi) - 5 (4V2) • Vi - 5 Vi • a|v2 + 4(/f2 - ATi) 

proving the first part. The second part follows from the definition ofK'. □ 

Proof of Proposition 16. Let y : [0, 6] — > M be the curvature line that is canonically associated with 
e, parameterized by arc-length and passing through p - 7(0). By definition, Cd is bounded above by 
the maximum distance from y to the tangent line passing trough y'(0). Since y is parameterized by 
arc-length, we obtain 

Cd ^ ' (26) 



where TCy denotes the maximum curvature of y as a space curve. Decomposing the curvature vector of 



y into its normal and geodesic components, and denoting by 7C" and the respective maxima of the 



norms of these components. Lemma 17 yields 



n<y<n<; + n<'y<n<+ max ^ , (27) 
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since TC" < 7C by definition of 7C. Consequently, we seek a lower bound for \6ii{-y{s))\. To do so, first 
observe that by definition of 7C' the derivatives of ki and K2 are bounded by 7C'. Hence, the function 6^ 
is Lipschitz with constant 27C', i.e., 

Mp)-6,(q)\<2'K'dM{p,q), 

for every point q e M. We now distinguish two cases: (i) |5a:(p)I < 4"7C"'e and (ii) ^ 4'K'e. In the 

first case, we immediately obtain 

Mp)ed\ < , 

which already proves the claim of the lemma. In the second case, we observe that for all s e [0, e] , we 
have 

1. 



Plugging this into (27) gives 



Together with (26), this yields 



\S.(7{s))\>\6,{p)\-2^'e>-\6,{p)\ 



\6ApM < \S,(P)\ I y + 1 £^ < {■K- + 7f ')<^ , 



completing the proof. 



3.4 Combining the strings 

We are now in the position to prove Theorem 2. 

To this end, assume that the edge e in our local polyhedral approximation is associated with the 
principal direction given by vi, and let this be the canonical direction if p is not umbilical. Furthermore, 
let - 2 sin y ||e|| be the discrete integrated curvature introduced in Section 2.1. Then propositions 12 
and 16 imply that there exists a constant C = C{'K, ,p) such that 



\K - 2A^K2\ < Ce^ 



(28) 



Assume additionally that e is chosen such that it maximizes the circumcentric area Ag. (There are exactly 
two choices for non-umbilical vertices.) Proposition 3 shows that this choice leads to lower and upper 
bounds for Ae by (1/C)e^ and Ce^, respectively, with C = C{p). Therefore, we can divide (28) by 2Ae 
to obtain 



2A, 



< Ce 



with C - C('K,'K',p). 

Finally, in order to prove our error estimate for the two other edge-based integrated discrete cur- 
vatures, ke = 2 tan yllell and = 6e\\e\\, we infer from Lemma 8 and a simple application of Taylor's 
theorem that there exists a constant C = C^K, p) such that 



dp - 2 sin 



< Ce^ and 



9e - 2 tan 



which can be used to obtain the requisite bound (28) for the other two discrete curvature definitions as 
well. This completes the proof of Theorem 2. 
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